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ABSTRACT 

Secondary CMB anisotropies and polarization provide a laboratory to study structure formation in 
the reionized epoch. We consider the kinetic Sunyaev-Zel'dovich effect from mildly nonlinear large-scale 
structure and show that it is a natural extension of the perturbative Vishniac effect. If the gas traces the 
dark matter to overdensities of order 10, as expected from simulations, this effect is at least comparable to 
the Vishniac effect at arcminute scales. On smaller scales, it may be used to study the thermal history- 
dependent clustering of the gas. Polarization is generated through Thomson scattering of primordial 
quadrupole anisotropies, kinetic (second order Doppler) quadrupole anisotropies and intrinsic scattering 
quadrupole anisotropies. Small scale polarization results from the density and ionization modulation 
of these sources. These effects generically produce comparable E and imparity polarization, but of 
negligible amplitude (10~ 3 — 1CP 2 fiK) in adiabatic CDM models. However, the primordial and kinetic 
quadrupoles are observationally comparable today so that a null detection of £>-polarization would set 
constraints on the evolution and coherence of the velocity field. Conversely, a detection of a cosmological 
B-polarization even at large angles does not necessarily imply the presence of gravity waves or vorticity. 
For these calculations, we develop an all-sky generalization of the Limber equation that allows for an 
arbitrary local angular dependence of the source for both scalar and symmetric trace-free tensor fields 
on the sky. 



1. INTRODUCTION 



With the rapid improvement in the sensitivity of cos- 
mic microwave background (CMB) experiments in recent 
years, it becomes increasingly important to address small 
contributions to the anisotropy and polarization at the fiK 



level and below. While the primary anisotropics and po- 



lariza tion from the epoch of recombination are thought 

to be understood at this level theoretically, secondary 
anisotropies from reionization are not. In this work, we re- 
visit the generation of these secondary anisotropies, uncov- 
ering a potentially important contribution from the non- 
linear regime, and explicitly calculate the small secondary 
polarization signal. We consider only those contributions 
from reionization that are true, frequency- independent, 
temperature distortions. Isolating these contributions ob- 
servationally from potentially larger but frequency-depen- 
dent foregrounds will be a great challenge b ut one that 
lies bey ond t he scope of this paper (see e.g 



Gispert 1999; Tegmark et al. 1999) 



Bouchct & 



Even the epoch of reionization itself cannot be accu- 
rately predicted in the context of an otherwise precisely 
defined theory of structure formation. The efficiency with 
which ionizing photons are produced and escape from the 
first baryonic objects in the universe will likely remain un- 
certain even with s ubstantial a dvances in numerical sim- 
ulations (see e.g. Abel et al. 1998). In the adiabatic 



cold dark matter (CDM) class of models, reionization is 
expected to occur in the range 8 <; zi <J 20 but with 



large uncertainties (see e.g. Haiman & Knox 1999 for a re- 
cent review). Observationally, reionization must be essen- 
tially complete by z ~ 5 due to the absence of the Gunn- 



Pcters on effect in quasar absorption spectra (Gunn & Pe 



terson 1965). For models with a roughly scale-invariant 



spectrum of initial fluctuations, detections of degree-scale 
anisotropies imply optically thin conditions or Zi <J 50 



( Griffiths et al. 1999). The large angle polarization of 
the CMB can in principle yield precise constraints on the 
reionization epoch but has not been detected to date. 
Even if reionization took place around z = 5, sec- 



ondary anisotropies from the Vishniac (1987) effect should 
produc e fiK anisotropies in the arcminute regime (Hu & 
Whitej 1996). The reason is that this is a second order 
effect whose amplitude per logarithmic interval in (1 + z) 
is constant, i.e. the contributions from below z = 5 are 
similar in magnitude to those in the range 5 < z <J 30. In 
the optically thin adiabatic CDM models, at least ~ 1/2 
of the Vishniac effect comes from redshifts z <; 5. Given 
these rather low redshifts of formation, it is interesting to 
consider whether nonlinear effects can further enhance the 
anisotropy. After developing general techniques for calcu- 
lating secondary anisotropies and polarization in §0, we 
consider contributions from the kinetic Sunyaev-Zel'dovich 
effect from large scale structure in the mildly nonlinear 
regime in §0, and show that it is the natural nonlinear ex- 
tension of the Vishniac effect. On small scales, both arise 
from the density modulation of the Doppler effect from 
large-scale bulk flows. 

We then turn to secondary polarization. Back-of-thc- 
envelope estimates in these optically thin conditions im- 
mediately place the secondary polarization signal orders 
of magnitude below the secondary anisotropies, except 
for the well-studied large-angle polarization from reion- 
ization. Nevertheless given the great potential of preci- 
sion polarization meas uremen ts for studying gravity waves 



(Kamionkowski et al. 1997; Zaldarriaga fc Seljak 1997), 



an explicit calculation is of interest. In we treat the 
polarization that results from second-order Doppler shifts 
due to bulk flows and separate the contributions in the 
two parity modes. This calculation may also be of inter- 
est to studies of CMB polarization in galaxy clusters. In 
§||, we consider the analogue of the Vishniac effect for po- 



1 



2 



larization: the density modulation of the kinetic and pri- 
mordial polarization sources as well as the generation of 
polarization through rescattering of Vishniac temperature 
anisotropies. Finally in §|| we show that the same rela- 
tionship between anisotropies and polarization for density- 
modulated signals holds for ionization-fraction modulated 
signals, as relevant for the brief epoch of inhomogeneous 
ionization expected at the onset of reionization. 

The basis for all of these calculatio ns is pr esented in 
the Appendix. There we generalize the Limber (1954) ap- 
proximation for sources with arbitrary local angular de- 
pendence and fields on the sky that are either scalar or 
tensor in nature. These formulae also extend the flat sky 
Limber approach to the full sky and may be useful in other 
cosmological studies. 

2. GENERAL CONSIDERATIONS 

We review the relevant properties and parameter s of the 
adiabatic CD M scenario for structure formation in §OI. In 
§[2.2| and §2.3, we discuss the general principles involved in 



calculating secondary CMB anisotropies and polarization 
based on formalism developed in the Appendix. 

2.1. Adiabatic CDM Model 
We work in the context of the adiabatic cold dark mat- 



ter (QDM) family of models where structure forms through 
the gravitational instability of ihe CDM in a background 

Friedmann-Robertson- Walker metric. In units of the crit- 
ical density 3H$/8ttG, where H = 100ft, km s _1 Mpc -1 is 
the Hubble parameter today, the contribution of each com- 
ponent is denoted f2j, i = c for the CDM, b for the baryons, 
A for the cosmological constant. It is convenient to define 
the auxiliary quantities Q m = fl c +Qi, and = 1— y"V fli, 
which represent the matter density and the contribution of 
spatial curvature to the expansion rate respectively. The 
expansion rate then becomes 



h 2 = Hi [n m (i + zf + n K (i + zf + n A ] 



(i) 



Throughout this paper, the notation A|. will always repre- 
sent the logarithmic power spectrum of the field "S" and 
should be assumed to be time-dependent even where that 
argument is suppressed. The rms of the field is defined as 



dk 

T 



-Ai 



(5) 



The cosmological Poisson equation relates the power 
spectra of the potential and density perturbations S 



Ho 
k 



ft 2 m (l + z) 2 A 2 , (6) 



and gives the relationship between their relative normal- 
ization 



(7) 



Here 5h is the amplitude of present-day density fluctua- 
tions at the H ubble scale; we adopt the COBE normal- 
ization for 5h ( Bunn fc White 1997). G(z) is the growth 
rate of linear density perturbations S(z) = G(z)S(0). Ex- 
pressions and approximations for this grow th fu nction can 
be found in the literature (e.g. Peebles! 1980; |Caroll ct 



al. 1992) and when the energy density in radiation can be 
neglected 

ow K m j~ 3 . (8) 

For the matter dominated regime where H oc (1 + z) 3 / 2 , 
G oc (1 + z)^ 1 and A = const. Likewise, the continuity 
equation relates the density and velocity power spectra 



a _ / GHq \ 2 
\ G k J 



(9) 



Convenient measures of distance and time include the 
conformal distance (or lookback time) from the observer 
at redshift z = in units of the Hubble distance today 
H^ 1 = 2997.9/i" 1 Mpc, 



D(z) 







dz . 



/o H(z') 

and the analogous angular diameter distance 
D A = n K 1/2 sinh(0^ 2 D) . 



(2) 



(3) 



Note that as il K -> 0, D A -> D. 

The adiabatic CDM model possesses a power spectrum 
of fluctuations in the gravitational potential $ 



a 2 _ p 



<i» 



.4 



k 

Ho 



T 2 {k) . 



(4) 



where the transfer function T(k) = 1 for scales much larger 
than the horizon at matter-radiation equa lity. We employ 
the CMBFast code flSeljak fc Zaldarriaga| 1996) to deter- 
mine T(k) at intermediate scales and extend it t o small 
scales using the fitting formulae of Eiscnstcin fc Hu| (1999). 



For this type of fluctuation spectra and growth rate, 
reionization of the universe is expected to occur rather 
late Zi <; 50 such that the reionized media is optically 
thin to Thomson scattering of CMB photons r < 1. The 
probability of last scattering within dD of D (the visibility 
function) is 



g = ie- T = Xt h (1 + z) 2 e- T . 



(10) 



Here dots represent derivatives with respect to D, X is the 
ionization fraction, and 



t h = 0.0691(1 -Y p )ft b h, 



(11) 



is the optical depth to Thomson scattering to the Hubble 
distance today, assuming full hydrogen ionization, where 
Y p is the primordial helium fraction. Note that the ioniza- 
tion fraction can exceed unity: X = (1 — 3Y p /4)/(l — Y p ) 
for singly ionized helium, X = (l — Y p /2)/(l—Y p ) for fully 
ionized helium. We assume that X = 1 in the reionized 
epoch such that 



■ 3Q m + y/1 + n m z(n m z + 3ft m - 2)] , (12) 
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Fig. 1. — Temperature anisotropies for the fiducial ACDM model 
with t = 0.1 (zj = 13). The secondary density (S^) modulated signal 
has been calculated under the assumption that the gas traces the 
dark matter. The ionization-modulated signal assumes patches of 5 
Mpc comoving size and duration of patchiness <5z;/(l + Zi) = 0.25. 
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Fig. 2. — Polarization for the fiducial ACDM model with r = 0.1 
separated into E (solid lines) and B (dashed lines) contributions. 
Secondary anisotropies from the primordial quadrupole are labeled 
(Prim. Q): (upper) homogeneous scattering; (lower) density (St,) 
and ionization (X) modulated scattering following Fig. M. For the 
kinematic quadrupole, the homogeneous and density modulated sig- 
nals are shown; the ionization modulated and intrinsic quadrupole 
signals falls below this range. Note that the B-parity polarization 
induced by gravitational lensing is much larger than any of these 
secondary B signals. 



for J7a 



and 

= 2 TH_ 
3 f^m 



y\ - n m + n m (i + zf - 1] 



(13) 



for tt K = 0. 

Although we maintain generality in all derivations, we 
illustrate our results with a ACDM model. The parameters 
for this model are tt c — 0.30, fi& = 0.05, £1a = 0.65, 
h = 0.65, Y p = 0.24, n = 1, and S H = 4.2 x 10~ 5 . This 
model has mass fluctuations on the 8h Mpc -1 scale in 
accord with the abundance of galaxy clusters erg = 0.86. 
A reasonable value here is important since many of the 



effects discussed below are nonlinearly dependent on the 
amplitude of mass fluctuations at or below this scale. We 
consider reionization redshifts in the range 5 < Z\ <J 30 or 
0.025 < t < 0.3. For reference, the primary anisotropies 
and polarization for a model with r = 0.1 [%i = 13) is 
shown in Figs. [I] and |^ respectively 

2.2. Secondary Anisotropies 

The motion of the scatterers imprints a temperature 
fluctuation on the CMB through the Doppler effect 



6(n) = J dDg(D)h-v b {*-), 



(14) 



where n is the direction on the sky, and Vf, is the velocity 
field of the baryons and is evaluated along the line of sight 
x = Dn. 

Let us first consider the case where the visibility g is 
dependent on time only. The source n • v& is a field that 
depends both on direction and on spatial position. It is 
useful to decompose the Fourier transform of the velocity 
field into a scalar and two vector (vorticity) modes 

v b (k)=-^)e 3 + E 4^-^, (15) 



m— ±1 



such that 



n.v 6 (k) = -»J^ «i m) n (m) (ft), (16) 



where e3 || k. As shown in the Appendix, the angular 
power spectrum resulting from the weighted projection of 
a source with a local dipolar angular dependence is (see 
equation | Al£ |), 



HH 



7T" 

1* 



dD g Da 



E 

m=±l 



A2(m) 



(17) 



for I 3> 1. The power per logarithmic interval in the ve- 
locity field A«f m ^ is evaluated at the wavenumber that 
projects onto the angular scale I at D, 



k = H Q 



D A 



(18) 



The implicit assumption here is that the visibility- weighted 
source g(D)vb(k, D) is slowly- varying across a wavelength 
of the perturbation. 

We will often use the power per logarithmic interval in 
[iK 2 or nK 2 = (10- 3 /iK) 2 



AT/ 



£(£ + V 

2vr 



CMB i 



(19) 



where T C mb = 2.728 x 10 6 /iK. The rule of thumb for AT 2 
is that it is of order the logarithmic power spectrum of the 
source (at cosmological distances Da ~ 1) divided by the 
multipole I. This factor of i comes from the loss of modes 
parallel to the line of sight due to crest-trough cancellation 
of the contributions. 



4 



Notice however that the m = potential flow compo- 
nent drops out of the final expression and violates this 
rule. In a potential flow waves perpendicular to the line 
of sight lack a velocity component parallel to the line of 
sight; consequ ently t here is no Dopplcr e ffect to leading 
order ( |Kaiser| 1984). pstriker fc Vishniac| (1986) pointed 
out that the same is not true for vortical flows since waves 
that run perpendicular to the line of sight have velocities 
parallel to the line of sight. Since flows in the linear regime 
are potential, the leading order Doppler effect is nonlinear 
in the perturbations at small scale. 

It is not necessary for the flows themselves to possess 
vorticity. Since it is the visibility-weighted velocity field 
that is the real source, spatial modulations in the vis- 
ibility create an effective velocity field 



v b (l + Sg/g) = v fc + ■ 



(20) 



which can be used in place of v b in equation ( |l7| ) . Spatial 
variations in the free electron density modulate the visi- 
bility and themselves may be caused by fluctuations in the 
net baryon density (see §0) or the ionization fraction (see 

§D- 

2.3. Secondary Polarization 

Thomson scattering of radiation with a quadrupolc 
anisotropy 

Q{m) ( x ) = - / dQ 2 A. J 9(x, n) (21) 



■in 

generates linear polarization in the CMB. In terms of the 
Stokes q and u parameters 

[ 9 ±z M ](n>=^ I dDg £ Q^ ±2 Y 2 m (n), (22) 



m=-2 



where Q^ m > are the 5 quadrupole moments of the tem- 
perature fl uctuation 1 and + 2YJ 71 are the spin-2 spherical 
harmonics (Goldberg ct al. 1967). 



Some subtleties arises when defining the Fourier trans- 
form of the quadrupole source since the orientation of the 
coordinate system enters into its definition. It is conve- 
nient to chose the reference frame so that 63 |j k, and 
this choice will be understood unless otherwise specified. 
The consequence is an ambiguity in combining contribu- 
tions from different fc-modes but these can be removed 
by defining t he coordinate indepen dent c omponents of the 
polarization ( Zaldarriaga & Seljak 1997; Kamionkowski et 
g 1997), 

[g± *«](*>=£ £ [E^±iB^] ±2 Y 2 m (n). (23) 

l m=-2 

Since polarization is a spin-2 field and its source has a 
local quadrupole ang ular d ependence, its power spectrum 
is given by equation ( A18| ) with s — ±2 and j — 2. The 
power spectra of the E and B parity states are then 



C EE -~- 



3tt 2 



lOf 3 



dDg 2 D A -A 



2(0) 

Q 



8 ^ 

m=±2 



1 Q( m ) = e^'/v^ in the notation of |Hu fc White) (1997), see 
their eqn. [77]. 




BB_ 



3tt 2 

"km 3 



[ dDg 2 D A \ £ A 

J 1 m=±l 



2 (m) 



(24) 



again evaluated with the projection equation ( |18| ) and a 
slowly-varying source. We define the logarithmic power 
spectrum of the polarization in the same way as for the 
anisotropies (see equation (l9)). It is likewise down by a 
factor of I from the power spectrum of the source except 
when the angular dependence of the source reduces it fur- 
ther as in the case of B-parity polarization from m = 0, ±2 
and _B-parity polariz ation fr om m = ±1. 

As point ed out by | Kaisci (1984) and examined quanti- 
tatively by Efstathiou| (1988), the polarization, even more 
than the temperature anisotropy, is suppressed by can- 
cellation at small scales. The reason is simply that the 
source of the polarization is the temperature anisotropy 
(quadrupole moment) and thus both the source and its 
contributions along the line of sight are suppressed. 

The dominant source at small angles is again the visi- 
bility modulation of the large scale quadrupole anisotropy 



Q 



(m) 



(1 + Sg/g) = Q {m) + 



)( m ) 



(25) 



We will consider quadrupole sources in §[|, the ef- 

fect of density modulation in §|5|, and that of ionization 
modulation in 

3. DENSITY-MODULATED ANISOTROPIES 

In this section we consider the temperature anisotropies 
that arise from the modulation of the Doppler effect due 
to density inhomogeneities. We begin with a review o f the 
second-order effect from linear density fluctuations (§ |3.l| ) 
and then generalize the calculation to include small-scale 
nonlinearities in the density field (§3.2). 



3.1. Linear Fluctuations: Vishniac Effect 

Spatial variations in the opacity due to density pertur- 
bations in the baryons modulates the visibility: 5g/g = 5^. 
When 5h is in the linear regime, the result is called the 
Vishniaq (1987) effect. In this section, we rederive the 
Vishniac effect in a manner that will make its nonlinear 
generalization obvious. 

The Fourier transform of v g is a convolution of the linear 
velocity and density fields, 



v s( k )=y ^p-vb(ki),5 6 (fc 2 ) , 



where here and througout 



k 2 = k - ki . 



(26) 



(27) 



The two vortical components to the velocity field are given 
by the projections (see equation Q) 



62 ± ie% 




(28) 



%(fci)5 b (fc 2 )F 1 ±1 *(k 1 ), 



in the basis where e3 || k. We have assumed that the un- 
derlying velocity field is a potential flow v b (k) = — w b (fc)k. 
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The power spectrum then becomes 

u3 



100 



A 2( 1) = l|*l^fe_ |yil(ki)|2[A + B]) 



(29) 



and likewise for the m = — 1 component. The two con- 
tributions are from the velocity-velocity, density-density 
power spectra 

A = Al b (k 1 )Al(k 2 ), (30) 
and the velocity-density cross correlation power spectra 

B = -^A 2 VhS M&sM, (31) 



where we have used the relation 

y 1 1 (k 2 ) = -|ir 1 1 (k 1 ). 



(32) 



Substituting equation (29) and (|17|) after reexpressing the 
velocity power spectrum with the density power spectrum 
using equation (O), we obtain 



dDD\lg^) Ajjy, (33) 



where we have used the fact that m = ±1 contributes 
equally. The mode coupling integral is 



Iv= 



vlvl 



Ai(fc) AUk) 



(34) 



Here and throughout 



VT&z/k = V 1 - 2/Ltyi 



(35) 



The form of the mode coupling integral in equation (34) 
is 1/2 that of Vishniacj (1987) equatio n (2.13), which rc- 
solves the appa rent discrepancy with Dodclson fc Juba^ 
(1995) raised by Jaffc fc Kamionkowski (1998). As pointed 
out by Dodelson & Jubas (1995), this form is easier to eval- 
uate numerically Additionally, for our purposes, it better 
brings out the small-scale limit. 

For wavelengths that are much smaller than the coher- 
ence scale of the velocity field, these expressions simplify 
considerably. The density- velocity cross correlation (31) 
vanishes and the remaining term ( |30| ) can be evaluated 
under the approximation k 2 = |k — ki| w k, 



aL (±1) = IaIv 



Si, "rms ' 



(36) 



where we have used the orthogonality of spherical harmon- 



ics / <MY™*Yf = S tti ,S m , m > and 



dk . n 

— A 

k Vb 



(37) 



A 2 / A 2 (lin) ; 




100 



(h Mpc 



Fig. 3. — Velocity-density decoupling and density nonlinearities. 
For the fiducial ACDM model the velocity and density fields decor- 
relate in the Vishniac calculation (Iy ec /Iv ~ 1) before the onset of 
nonlincarity in the density field, especially at high redshift where 
most of the scattering occurs. 



Note that v 2 a is still a function of D. 

Equation ( pq ) has a simple interpretation. At small 
scales, the effect arises from a density modulation of a 
uniform bulk velocity from larger scales whose rms in each 
component is 1 /3 of the total. 

Substituting the velocity power spectrum relation (36) 
into equation (j2^) yields 



2tt 2 



2 „,2 

rms ' 



(38) 



This equation may alternately be derived from equation 
( |33| ) in the limit that the integral gets its contribution 
from yi <§; 1, 

(dcc) _4 f 00 dy 1 A%{ky u rj) 



1 T/ — "T 



3 jo 

4 (G_k_ 
"3 \GH~o 



Vi A^(fc,r?) 



k2 

2 2 
V 



(39) 



In Fig. ||, we show that for the ACDM model, the density 
and velocity contributions decouple in this manner beyond 
k ~ 0.2/iMpc -1 . The approximation ( (39|) represents an 
efficient way of evaluating the otherwise computationally 
expensive the mode coupling integral in this regime. 

3.2. Nonlinear Perturbations: Kinetic SZ Effect 

Nonlinear density fluctuations caught up in a bulk 
flow from larger scales gives rise to the kinetic Sunyaev- 
Zel'dovich effect from large-scale structure. It can alter- 
nately be thought of as the nonlinear extension of the Vish- 
niac effect. The reason is that in adiabatic CDM cosmolo- 
gies nonlinearities only affect the density field below the 
coherence scale of the bulk velocity. Recall that in this 
limit the Vishniac effect arises from density perturbations 
caught in a large-scale flow to which it is uncorrelated. 

Figure H illustrates this fact for the fiducial ACDM. No- 
tice that even at z = 0, the dark matter density field only 
deviates substantially from the linear approximation after 
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Nonlinear Power Spectra 



Fig. X 

der the 
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Nonlinear v« linear power spectra for the CDM un- 



Peacock & Dodds (1996) scaling approximat ion and for the 
baryons under the additional prefiltcring ansatz of Gncdin & Hui 
(1998). Note that the dilation of scales brings the filtering scale 
deep into the nonlinear regime. 
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Fig. 5. — Maximal nonlinear enhancement of the Vishniac effect. 
Under the assumption that the gas density traces the dark matter 
density into the deeply nonlinear regime the Vishniac effect is sig- 
nificantly enhanced by nonlinearities at i > 1000 especially in the 
late reionization scenarios. 



equation ( J39| ) becomes a good description of the Vishniac 
effect. If we then replace the linear density power spectrum 
with its nonlinear analogue but leave the contribution from 
the velocity power spectrum the same, we obtain 



ee 



7T 



dDD\ 



G 
'G 



2(li: 



(40) 



with [he mode coupling iuLegral I v using equation (34f 

evaluated under linear theory. This expression has the 
nice feature that it expresses the total effect: it includes 
both the Vishniac effect and the kinetic SZ effect from 
nonlinear structures. 

The underlying assumption is that the density fluctu- 
ations in the nonlinear regime are uncorrelated with the 
bulk velocity field. Nonlinear evolution in the density field 
will correlate velocity and density modes if they are both 
in the nonlinear regime. However, the bulk flow in adia- 
batic , CDM models arises mainly from the linear regime. 



In thd A CDM monVl ; half the contributions to vf mti comes 



from scales k < O.OTfeMpc -1 and the fluctua tions go non- 
linear at k ~ 0.2/iMpc" 1 . Scoccimarro et al. (1999) found 
that for the 4-point statistics, nonlinear modes are highly 
correlated but that the correlation drops rapidly as one 
pair of the modes enters the linear regime. Decorrelation 
is even more effective here since the relevant velocity and 
density modes are oriented perpendicular to each other. 

Once the nonlinear power spectrum of the baryonic gas 
is known, the resultant anisotropics can be calculated with 
equation (f40|). This unfortunately requires hydrodynamic 
simulations i n general a nd ab initio attempts in the lit- 
erature (e.g. Pcrsi ct al] 1995) to simulate the scattering 
effects have not had the dynamic range to treat the full 
problem. Let us instead break the problem into two: cal- 
culate the effect under the assumption that the gas tracks 
the dark matter and then estimate where this approxima- 
tion may break down. 

3.2.1. Maximal Effect 



Let us first consider the effect under the simple assump- 
tion that the gas traces the dark matter to place an upper 
limit on the magnitude of the effect. The nonlinear dark 
matter power spectrum has been well-studied with N-body 
simul ations for the fu ll range of adiabatic CDM cosmolo- 
gies. Hamilton et al.| (1991) suggested a useful scaling re- 
lation for the correlation function in the nonlinea r regime 
that was generalized to the power spectrum by Peacock 
& Dodds (1994). The underlying idea is that nonlinear 



fluctuations on a scale k arise from linear fluctuations on 
a larger scale 



AUk)]- 1/3 k- 



(41) 



so that there is a functional relation between the nonlinear 
and linear power spectra at these two scales 



Al a (fe) 



/ NL [Af in) (fc ni 



)] 



(42) 



which may be fit to simulations. We take the Peacock & 
Doddsj (1996; eqn. [21]-[27]) form for / NL . The resulting 



power spectrum for the ACDM model are shown in Fig. kt[ 
Under this assumption, the kinetic SZ effect provides a 
significant enhancement of the Vishniac effect. As shown 
in Fig. H the relative contribution is the largest for late 
reionization scenarios (low z,-) because anisotropies then 
arise from structure that is well-developed. However one 
must keep in mind that this is essentially an upper limit 
to the contributions since the gas pressure will smooth out 
the gas density below the Jeans scale. 

3.2.2. Pressure Cutoff 



Gncdin & Hui (1998) examined simple schemes to ap- 
proximate the effect of gas pressure. One such scheme that 
has fractional errors on the 10% level for overdensities <; 10 
is to filter the density perturbations in the linear regime as 
A$ = f£(k/kp)A'fi a and treat the system as collisionless 
baryonic particles. Their best fit is obtained with the filter 



fb 



(l+4x 2 ) 1 /4 J 



(43) 
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Fig. 6. — Prefiltering the dark matter spectrum. Pmfiltering the 
power spectrum to account for gas pressure as in Fig. rN cuts off the 
anisotropies at £ ~ 10 4 for the ACDM model. 
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Fig. 7. — Postfiltcring the dark matter spectrum. Filtering the 
CDM power spectrum with equation QMJ) introduces a smooth cut 
off at small angles and shows that nonhnearities are important as 
long as the gas traces the dark matter out to Ap ~ 10. 



1/2 

and Gncdin (1998) suggests kp = 34Qm /iMpc _i as a rea- 
sonable choice for the thermal history dependent filtering 
scale. 

The results of applying the nonlinear scaling relation 
of equation (B2| to the filtered linear power spectrum are 
shown in Fig - Hi Notice that the dilation in wavenumber 
takes the filtering scale deep into the nonlinear regime. 
The calculations should not be trusted beyond this scale 
since the nonlinear scaling breaks down for spectra with 
a sharp cut off; indeed we have taken the local slope of 
the unfiltcred spectr um rather than the filtered spectrum 
when evaluating the Peacock fc D_odds| (1996) formulae. 
Still the result displayed in Fig. || imply that the gas- 
traces-CDM assumption is reasonable in the arcminute 
regime £ <; 10 4 . 

To give a more model-independent quantification of this 
effect, let us also consider a baryonic power spectrum given 

by 




0.0001 



(h Mpc - 



FlG. 8. — The power spectrum of the present-day quadrupole 
from the projection of primordial temperature inhomogeneities and 
the kinematic Doppler effect in the ACDM model. The primordial 
quadrupole has been calculated numerically with a Boltzmann code. 



The results for several values of A^ are given in Fig. ^ and 
show that only if the baryons fail to trace the dark matter 
out to Ap ~ 10 can nonlinear effects be ignored. 

Our analysis shows that the kinetic SZ effect should 
be an important contributor to small- angle anisotropies. 
On arcminute scales and above, the assumption that the 
gas traces the dark matter is reasonable given our cur- 
rent understanding of the thermal history of adiabatic 
CDM models. The amplitude of the effect at the sub- 
arcminute scales will depend on the details of thermal his- 
tory throughits effect on the clustering of the gas. Ob- 
servations in this regime may open a new window on the 
physics of this regime. Isolation of this effect from fore- 
grounds (e.g. radio and infrared point sources) and other 
secondary anisotropies (e.g. the inhomogeneous reioniza- 
tion signal considered in §o) will be the main obstacle to 
overcome. 

4. POLARIZATION FROM HOMOGENEOUS SCATTERING 

In this section we consider the origin of quadrupole 
anisotropies and the linear polarization they generate 
through Thomson scattering in a homogeneous medium. 
The quadrupole anisotropies fall into three broad classes: 
the primordial quadrupole from the projection of Sachs- 
Wolfe temperatu re a nisotropies, the intrinsic quadrupole 
from scattering (§4.1), and the ki nem atic quadrupole from 
the second order Doppler effect ( §|4.2| ). 



4.1. Primordial and Intrinsic Quadrupoles 

A large-scale temperature inhomogeneity at the recom- 
bination epoch (z ~ 10 3 ) from the Sachs- Wolfe gravita- 
tional redshift effect O = — $/3 becomes a quadrupole 
anisotropy at low redshifts by simple projection. In other 
words, the spherical harmonic decomposition of the plane 
wave fluctuation, 



4> 



exp(ik ■ x) 



6 Hq 



(45) 



Ai(fc) = Al(fc)exp(-Al/A£) 



(44) 
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yields a spectrum of anisotropies and in particular a 
quadrupolc anisotropy given by 



(46) 



Here D* = D ICC — D and we have assumed that spatial 
curvature can be neglected at reionization. Curvature cor- 
rections may be included by replacing the spherical bessel 
function fa with the hyperspherical bessel function $2 (see 
Appendix) . 

The logarithmic power spectrum of the quadrupole can 
therefore be expressed in terms of that of the potential as 



A^i 2 2 (^*)A 2 



k 



(1 + -) 2 G 2 ^1(^), (47) 

where the second line employs equation (Q) and again as- 
sumes matter domination at D. We show the results of 
numerically calculating the quadru polc power spe ctrum 
from a cosmological Boltzmann code (White & Scott 1996) 
in Fig. ^. Note that the spectrum peaks at large scales, 
corresponding to the horizon at the reionization epoch 
k ~ Ho/D*. On smaller scales, the projection ( fig ) car- 
ries the power to higher angular moments. 

The rms is the integral over the power spectrum 



2 _[ * A 2(0) 



4^(1 + zfG^KD.Ho) 1 -" I dxx n jl{x), 



=-^(1 + z) 2 G 2 ^( J D*ff ) 1 -'Tsw(«) , (48) 



where 



^ r[(3-n)/2]r[(3 + n)/2] 

r sw (n) - 3 ^ r[(4 _„ )/2]r[(9 _ n)/2] 



(49) 



For reference, note that Q 2 ms = 5Cf )e /47r. For our fiducial 
ACDM model 6 H = 4.2 x 1(T 5 , (1 + z)G = 1.24, and 
rsw(l) = 1) we obtain QrmsTcMB = 16/_tK. 

There are two effects that modify this result slightly. 
The first is that the universe may not be completely 
matter-dominated at the time of scattering. Additional 
anisotropies are then created as the gravitational poten- 
tial decays. At high redshifts, the effect of the radi- 
ation energy density pushes the quadrupole anis otropy 
up vi a th e "early" integrated Sachs- Wolfe effect ( Hu & 
Sugiyima 1995); correspondingly at low redshifts the ef- 
fect of the cosmological constant or curvature again pushes 
the qu adrupolc up via the "lat e" integrated Sachs- Wolfe 
effect ( Kofman fc Sta ro binskii 1985). These effects are 
clearly visible in Fig. |ic|, where we plot the evolution of 
the rms quadrupole from numerical calculations. 

The second effect is that scattering itself can produce 
quadrupole anisotropics through projection of the Doppler 
effect and also destroy the primordial quadrupole through 
isotropization. We call the net effect the contributions of 
the intrinsic quadrupole. This effect is suppressed both 
by the assumed optically thin conditions and by the fact 



that the Doppler effect is cancelled except on scales near 
the horizon where the velocity itself is small. In Fig. [l0[ 
this effect is barely visible as a blip near the redshift of 
reionization. 

The hallmark of all these effects is that they create only 
m = type quadrupoles, because they come from [1 = 0, 
m = 0) gravitational potential perturbations and (1=1, 
m = 0) potential flows. They therefore generate only im- 
parity polarization. The polarization from these sources 
are not slowly varying across a wavelength of the horizon- 
sized perturbations and hence the approximation devel- 
oped in the Appendix cannot be used to calculate these 
effects. They are however automatically included in cos- 
mological Boltzmann codes. We show the resulting polar- 
ization in Fig. U for the ACDM model as calculated with 
CMBFast. As is well known, the secondary polarization 
from these effects are confined to the lowest £ since the 
quadrupole sources contribute mainly on horizon scales at 
last scattering. 

4.2. Kinematic Quadrupole 



As pointed out by Sunyacv fc Zel'dovich (1980), in the 
rest frame of the scatterers, an isotropic CMB gains a 
quadrupolc anisotropy from the quadratic Doppler effect 



6 = 



1 - n • Vb 



- 1 



n • v fa + (n • v 6 ) 2 - i«g , (50) 



that induces a polarization in the CMB by Thomson scat- 
tering. In real space, the quadrupole moment is simply 



O (0) (x) 



3^5 



(51) 



in the basis aligned with the velocity e 3 || vt(x); in Fourier 
space with a basis £3 || k, the expression is more involved 



n- v b (ki)n- v 6 (k 2 ) 

(27T) J 




Yi m *(ki 



(52) 



where = 2 = 0. The logarithmic power spectrum 
becomes 

C^^^^^^' (53) 
where the mode coupling integrals are given by 



Here the angular arguments are 



(54) 



a( ±2 ^y 2 (l- M 2 ) 2 , 
a(±^(l-2 Wl ) 2 (l-M 2 ), 



(55) 
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Fig. 9. — Power spectra of E- and B-parity polarization from the 
kinematic quadrupole in the ACDM model. 



and it is useful to define the auxilary quantity 
a (/) = |(l-M 2 )- 



(56) 



The polarization power spectrum then follows by insert- 
ing equation ( |53| ) in equation (^Jj). The result for the 
ACDM cosmology is shown in Fig. g 

The rms quadrupole follows simply from the relation 



£ « (m) = yl a 



(57) 



m=-2 



with a change of variables to k2 for the integration of the 
first term 



) 2 = 

! rms 



/kin 



dk 
~~k 



2 



n=-2 



2 (m) 



= ^(l-/kinX s (58) 



1 fdkAiAi 



r(/) 



For the ACDM model /kin = 0.16. This expression 
is smaller than the naive expectation implied by real 
space relation (51), ^{vf) /45 = 12w^ ms /45, because the 
quadrupole is oriented and not simply a scalar gaussian 
random field. The main effect can be understood as the 
contribution from the variance of v%, (vf\ — (yl) = 2vf ms , 
since its mean value comes from contributions that cancel 
out in orientation for the quadrupole. 

The kinematic quadrupole is not negligible in compar- 
ison with the primordial quadrupole for ACDM or in- 
deed any model that fits the observations. Large angle 
anisotropies are observed today at 1CP 5 while the veloc- 
ity field reaches 1CT 3 . As shown in Fig. [l0[ for ACDM it 
is approximately 1/4 the primordial one in amplitude to- 
day. Furthermore, equation ( |53| ) implies that it generates 
comparable power in all 5 m-modes of the quadrupole and 
hence comparable E and B parity polarization. One might 
naively assume that the kinematic effect generates as much 
£?-parity polarization as the primordial quadrupole gen- 
erates E. This would then be an obstacle for detecting 
gravity waves through the B-parity polarization. 
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Fig. 10. — The time evolution of the rms quadrupole in the ACDM 
model. Upper curves: the various linear theory effects (see text) 
for three different ionization histories from numerical calculations. 
Lower curves: the kinematic effect. 



There are however two reasons why the polarization is 
much smaller than implied by these arguments. The first is 
that the v^ ms declines as (l + z) -1 in the matter-dominated 
epoch whereas the primordial quadrupole remains con- 
stant. The second is that in the ACDM model the co- 
herence of the velocity field is over 2 orders of magnitude 
smaller than that of the primordial quadrupole and hence 
generates polarization at I ~ 400. Recall that the rule 
of thumb for projected sources is that the angular power 
spectrum is down by one factor of I from the spatial power 
spectrum. This produces another order of magnitude sup- 
pression of the polarization (see Fig. ^|). The .B-parity po- 
larization generated here is significantly smaller than that 
gener ated by gravitational lcnsing of the primary polariza- 
tion (Zaldarriaga & Seljak 1998). This effect is therefore 
unlikely to be detectable in ACDM models. 

It is worth emphasizing that the amplitude of this effect 
is highly model dependent and scales in power as vf ms X v 
where A„ is the coherence scale of the velocity field. Up- 
per limits on the B-polarization can therefore constrain 
the amplitude and coherence of the velocity field. For ex- 
ample consider a velocity field with an rms of ~ 700km s _1 
on the a tophat scale of \5 0h~ 1 Mpc is suggested by the 
data of |Laucr fc Postman (1994). If we model that field 

Ax 



by the ACDM velocity power spectrum but increase the 
coherence length by a factor of 10 and the rms velocity by 
1.55, we boost the polarization power by 1.55 4 x 10 = 58 
and angular scale by 10 (£ ~ 40). The power would be 
even larger in an open universe where the velocity am- 
plitude actually grows substantially with redshift before 
declining. It is also worth bearing in mind that effects 
such as these warn against blindly taking even a large- 
angle detection of S-polarization as a model-independent 
detection of gravity waves or vorticity. 

Finally, note that the correlation between primary and 
secondary temperature anisotropies and the kinematic po- 
larization is suppressed since it involves the three-point 
function of the density field and hence vanishes in the lin- 
ear regime. 
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5. DENSITY-MODULATED POLARIZATION 

Here we consider the polarization generated from the 
density-modulation of the quadrupole anisotropies con- 
sidered in the last section. Again these are separated 
by the natur e o f the quadrupole source: th e p rimordial 
quadrupole §^], the kinematic quadrupole §5.2, and the 
intrinsic quadrupole §|5.3| . 

5.1. Primordial Quadrupole 

For the modulation of the primordial quadrupole by 
density fluctuations, we have 



Qi m \k) 



Q (m) (ki)fc(*a), 



(59) 



Scalar perturbations in linear theory only generate m = 
quadrupoles in the ki-basis, i.e. Q( Q )(&i)Y 2 (/?, a) where 
P and a are the polar and azimuthal angles defining n in 
this basis. One can project this onto the e3 || k ba sis with 
the help of the angular addition relation (sec Bu 
1997 eq. [7]), 



White 



Y 2 °((3,a) 




(ki)KT(n) 



(60) 



such that 



d 3 ^ 
(2^)3 



)(0) 



(fci)Y2r(ki)<W , (61) 



Since the quadrupole source '(ki) peaks on the scale of 
the horizon whereas the density perturbations rise to small 
scales in adiabatic CDM models, the correlation decouples 



A 



2(m) 



a 2 n 2 



(62) 



where we have again used the orthogonality of spherical 
harmonics and Q 2 ms is taken from numerical calculations 
as in Fig. |l^. Inserting equation ([32]) into equation (pij), 
we obtain 



EE 



BB 



3tt 2 
50^ 3 



dDD A g 2 A 2 Q 2 , 



(63) 



As shown in Fig. ffll even for the gas-traces-CDM assump- 
tion the polarization generated by this effect in a ACDM 
model is small, much smaller than the polarization gener- 
ated by lensing (cf. Fig. |2|) . 

The cross correlation between multipolc moments of the 
Vishniac effect and the polarization Cf E vanishes when 
averaged over the sky due to the different angular depen- 
dence of the underlying dipole and quadrupole sources. 
Mathematically, this can be seen in the sources © and 
( pgj ) ; the orthogonality of the spherical harmonics guaran- 
tees that these terms will integrate to zero when summing 
over the directions of each A:-mode. 

Note however that there will be a strong correlation be- 
tween the amplitude of the polarization and the amplitude 
of the temperature fluctuations and may provide a means 
of pulling out the signal in models where it is somewhat 
stronger. 
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Fig. 11. — Density-modulated polarization in each parity type 
E = B for the ACDM model. Here we assume that the gas traces 
the dark matter. 



5.2. Kinematic Quadrupole 

The kinematic quadrupole is also modulated by the 
density field. Its effect is identical to the primordial 
quadrupole except equation (|5q ) for the kinetic quadrupole 
is used in equation (163) , 



/~iEE _ /-lilts 



BB 



4tt 2 
375£ 3 



J dDg 2 D A A 2 b (l-f kin )v. 



4 

rms > 



(64) 



where recall that /km 1S a small correction due to the direc- 
tionality of the quadrupole moments. The ratio integrands 
between the modulated Doppler and cosmic polarizations 
is 



45 



(1 - /kin)"^ 2 
Vrn 



(65) 



Since observationally i> rms ~ (1 + z) _1 ^ 2 10 -3 and Q rms ~ 
10 -5 , the effect from the primordial quadrupole is typically 
larger. We show the results for the fiducial ACDM model 
in Fig. | 

5.3. Intrinsic Quadrupole 

The quadrupole anisotropy generated by the Vishniac 
effect during reionization produces a linear polarization 
intrinsic to the Vishniac effect. This effect is doubly sup- 
pressed by cancellation: both the source of the quadrupole 
and its polarization contributions along the line of sight 
cancel when integrating along the line of sight. 

The quadrupole moment generated by the Vishniac ef- 
fect is given by equation (|A6j) 



and thus the power per logarithmic interval 

.2 



A^L|LW X ^(1 + ^A 



2 (m) 



(66) 



(67) 
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(a) Temperature 




Fig. 12. — Ionization-modulation effect. A patch size of 5Mpc and 
duration of patchiness of Szi/(1 + zi) = 0.25 in the ACDM model is 
assumed. 



The polarization power spectrum becomes 



„A 2 ( lin )A 2 T 



(68) 



with Cf E = since only m = ±1 contributes to the Vish- 
niac effect. The integrand in this equation is down by a 
factor of 



9?H 
AOOi 2 



D' A X 2 T^(l + zf = 



r r 2 



(69) 



compared with the temperature anisotropies and is thus 
highly suppressed at high I ~ 10 3 from an amplitude that 
is already reduced by the low optical depth t <J 1CP 1 in 
the reionized universe. Finally, there is no cross correlation 
between the temperature and intrinsic polarization since 
the polarization is purely of S-parity. 

6. PATCHY REIONIZATION 

Inhomogeneities in the ionization fraction create a mod- 
ulated Doppler anisotropies and polarization of the CMB 
just as inhomogeneities in the density. The techniques 
used in the calculation of density-modulated effects thus 
can be carried over to here with little modification. 

For the modulated doppler effect, one takes v g of equa- 
tion ( pp| ) to be VbSx where Sx is the fluctuation in the 
ionization fraction. Once the power spectrum of this quan- 
tity is known, the anisotropies follow directly from equa- 
tion @. 

Unfortunately, even the crude form of the power spec- 
trum is not known as it requires not only an understanding 
of the first baryonic objects in the unvierse but also a full 
three dimensional radiative transfer in an inhomogeneous 



medium. Thus even though we would expect the ionizing 
sources to be associated with the peaks in the density field, 
it may be that the radiation escapes into and ionizes the 
low den sity medium where th e recombination rates are the 
lowest (Miraldc-Escudc et al 1999). The techniques intro- 
duced here are therefore more useful for the inverse prob- 
lem. If such a signal is detected in the CMB, we will want 
to use these relations to invert it to uncover the power 
spectrum of the ionization fraction and hence learn about 
the manner in which the universe was ionized. 

Nevertheless it is us eful to have a co ncrete model to il- 
lustrate these effects. Gruzinov fc Hu (1998) introduced 



a toy model in which the universe is reionized in uncor- 
rected spherical patches of a characteristic size R in a 
rcdshift interval Az around Zi. A distribution in R, like 



the one considered by Aghanim et al. (1996), can of course 
be constructed by superimposing simple sources. It is im- 
portant to bear in mind that in a more realistic scenario 
there arc likely to be correlations between the ionized re- 



gions; 



Knox et al. 



(1998) work out the consequences of 
ionized regions tracing the peaks in the density field. In 
the simple uncorrelated model, the correl ation fuction of 
the ionization fraction can be modeled as (Gruzinov & Hx 
1998) 



(,5A'(x 1 ) < 5X(x 2 )H^ - X 2 ]e 



-[( Xl -x 2 ) 2 /2fl 2 ] 



(70) 



where X{D) now is the mean ionization. This form is 
chosen to have the right asymptotic behavior in time and 
scaling with R; other forms can be chosen which essentially 
correspond to a r edefinition o f the patch size R (e.g. the 
tophat spheres of Knox et al. 1998) 



With this ansatz, the logarithmic power spectrum be- 
comes 

X 2 A 2 Sx = ^(kRfe- k2R "/ 2 [X - X 2 } . (71) 

The angul ar p ower spectrum that results is simply that of 
equation ( p8| ) with the density power spectrum replaced 
with the ionization power spectrum 

(72) 



(2tt) 3 / 2 



3 / X 2 



2 [X-X 2 ]v 2 n 



where 9 = HqR/D 



Following Gruzinov & Hu (1998), we take a mean ion- 
ization fraction that grows linearly from zero at Di up to 
unity by Di + SDi 

D-D, 

Assuming reionization occurs promptly (SDi/Di -C 1), the 
other quantities may be taken out of the integral leaving 
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Likewise, the polarization follows from equation (|6 
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(76) 



with Qrms being cither the primordial quadrupole of equa- 
tion ( ^q) or the kinematic quadrupole of equation ( |58| ) . For 
adiabatic CDM models, the former typically dominates the 
latter even more so than for the density- modulated effect. 
The reason is that one no longer is weighted toward low 
redshifts and larger velocities by the growth of density 
perturbations. An example of the temperature and polar- 
ization signals for a specific choice of parameters is given 
in Fig. Ill 



7. DISCUSSION 

We have explored the role of mildly-nonlinear density 
fluctuations in generating secondary CMB anisotropics 
and a host of contributions to the secondary polarization. 
In adiabatic CDM models, the Doppler effect from scat- 
tering off nonlinear baryonic clumps in a large-scale bulk 
flow is a natural extension of the Vishniac effect. In the 
small scale limit, the Vishniac effect simplifies to the corre- 
sponding effect for linear density perturbations. For this 
class of models, the nonlinear contributions out to over- 
densities of ~ 10 are comparable to or greater than the 
linear contributions, depending on the redshift of reion- 
ization. If the gas traces the dark matter to much higher 
overdensities then upto an order of magnitude increase in 
the subarcminute anisotropy would result. Likewise small- 
scale patchiness in the ionization can greatly increase the 
anisotropy in this regime. Calculations of the dumpiness 
of the ionized gas are difficult and unreliable even with 
state-of-the-art numerical simulations. An observational 
study of subarcminute scale anisotropies therefore offers 
the opportunity to discover aspects of reionization that 
are intractable to theoretical analysis today. 

The secondary polarization from reionization is, as ex- 
pected, extremely small in the context of adiabatic CDM 
models for structure formation and unlikely to inhibit ex- 
traction of even subtle effects such as the -B-parity polar- 
ization from gravitational lensing or gravitational waves. 
It is worthwhile to note that these mechanisms generically 
predict comparable power in E and B parity polariza- 
tion and may be important outside of the adiabatic CDM 
model context. For example, in a model with a coherent 
1500km s _1 velocity field on the 100/i _1 Mpc scale the kine- 
matically induced quadrupole generates a .B-parity polar- 
ization at the 0.1/iK level for reasonable optical depths 
t ~ 0.1 — 0.3. Likewise, if the velocity field fell off less 
rapidly than (1 + z)" 1 / 2 out to redshifts of order 10 or 
more, then the effect might also be enhanced to this level. 
In these models, the density and ionization modulation of 
the scattering produce even larger signals in the arcminute 
regime that may be observable. Thus even a null detec- 
tion of .B-parity polarization can be interesting. The tools 
that we have introduced in the text and appendix are quite 
general and have applications beyond the adiabatic CDM 
model. 

Even in the adiabatic CDM class of models, signifi- 
cant uncertainties remain expecially in the contribution 
from inhomogeneities in the ionization and the small-scale 
behavior of the gas. These questions will ultimately be 
answered by the observations themselves as CMB experi- 



ments probe the arcminute regime with higher and higher 
sensitivity. 

APPENDIX 
GENERALIZING THE LIMBER EQUATION 



The Limber (1954) equation describes the two-point 
statistics of a field which is the two-dimensional projection 
on the sky of a three-dimensional source field whose statis- 
tical properties vary slowly along the line of sight. Kaiser 



(1992, 1998) expressed the result as a relation between the 



two and t hree dimensional power spectra of the fields. Hr 
fc White] (1996) generalized these flat-sky approximations 
in an all-sky approach for spatially flat cosmologies. 

Here we extend the techniques in three ways. We gener- 
alize them to open cosmologies where Fourier analysis on 
the three-dimensional field is invalid. We allow for an arbi- 
trary angular dependence in the source. Finally, we treat 
tensor fields on the sky in addition to the usual scalar 
fields. 

Total Angular Momentum Method 

A general field X that depends both on position x and 
direction n at x can be expanded in a complete set of 



modes denoted S G™ (see |Hu et alj 1998) 



X(x,n) 



d 3 c 



(2tt) 3 ^ 
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where the spin index s = 0, ±2 for scalar and tensor fields 
on the sky and (£,m) describe the multipolc moments of 
the local angular dependence. In flat space, these modes 
are simply the product of p lane waves and spin-spherical 
harmonics ( Goldberg et al. 1967) 



sGT = HY 



■in 



21+1 



.^(ftjexpfiq-x). (A2) 



Note that qY™ = Y™, the ordinary spherical harmonics. 
The angular power spectrum of the field is defined as 

cT-^J^^r. (A3. 

Let us suppose that the field on the sky X is generated 
by the line-of-sight integral of another positionally and 
directionally dependent field S evaluated at a x = Dh 

X(x,n)=J dDS(x,n) , (A4) 

=/<«>/ (^E^r- 

v 7 jm 

The normal modes can be rewritten in spherical coordi- 
nates as 



* G ? = EHJVM^+lJag sYl 



(A5) 



where we have used the spherical harmonic decomposition 
of a plane wave to combine the local and plane- wave angu- 
lar dependence terms into the total angular dependence as 
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seen by the observer at the origin. Here the o^™ 1 ] are lin- 
ear combinations of spherical Bessel functions with weights 
given by Clebsch-Gordan coefficients as we shall see. The 
moments of the field are then 



Xt } = [dD(2e + l)J2s( m) a ( ™ 



(m) ~ (m ] . (A6) 



For open geometries, the expressions in spherical coordi- 
nates take the same form except the radial functions & 
are linear combinati ons of the hyperspherical Bessel func- 
tions <f>^ (flu et al. 1998). They depend separately on 
radial distance 

(A7) 



X = ^k D - 



and wavenumber, 



(A8) 



which itself differs from the usual eigenvalue k of the Lapla- 
cian near the curvature scale 



q=Jk 2 - (|m| + l)Hfo K . 



(A9) 



Useful radial functions include 

..(0) . 




1(1+1) 
2{v 2 + 1) 



cschx . 



(2) _ 3 (i + 2)(P-l)£ 



(o)_l 

*02,r>> 



8 [v 2 +4)(i/ 2 - 


hi) 






/ i 


y (y 2 +4) (^ 2 - 


hi) 


\% £(£ + 1) 



csch 2 x §\ 



l > 



^1 (^+4x^+1) i^Jwr • (mo) 



for 



and 



^±22,1 



/3 + 


'2)(^ 2 ~ 


IK 


' 8 (y 2 


+ 4)(z/ 2 


+ 1) 


1 { " 


-!)(* + 


2) 


]j(u 2 


+ 4) (1/2 


+ 1) 


+$?'] 






/ i 


V 


+ 4) (l/ 2 


+ 1) 



csch 2 x$£ , 

cschx[(cothx ± iv)^ v t 

[$f' + 2(2cothx± 

+ (2coth 2 x + 1 - v 2 ± 4ii/cothx)$^] , (All) 
for s = 2. Note that primes are derivatives with respect 



to x and 



a 



Weak- Coupling Approximation 

The we ak-coupling approximation was introduce d in [Hu| 
fc White (1996) as a means of evaluating equation (|A6|) in 



the limit that the source sj"^ varies in distance D only on 
a much larger scale than the wavelength q^ 1 . In this case, 
the source may be taken out of the integral and evaluated 
at 



LH1 



~ 1/2 sinh- 1 ^ 2 ^). 



K 



(A12) 



The remaining integral over o^™"\ becomes 

(m) _H ^ T[(£+l)/2] (m) 

^V~r[(^ + 2)/2]^ x " (A13) 

where the correction due to spatial curvature is given by 
K v = [l+ (£+1 / 2) V 1/4 tanh 1 / 2 (g,). (A14) 



and the weights wjjj 1 ^ follow from the radial functions 



(A10) and (All) 



W {0) -\ 
1 1 oo —■ L 




csch 2 x , 



(0)_ 



02 ■ 



1 / (^2 + l) 

= 2V(^ 2 +4)' 



(A15) 



for the spin zero fields and 



wi°^l- 



w. 



(1)_ J 

±22~ r 



w. 



(2)_ J 
±22-, 



/3 (^ + 2)(£ 2 - 


l)i 


' 8(i/2+4)(^ 2 




/ V-W + 


2) 


V (f 2 +4)(i/ 2 


^ 1) 


x (cothx ± 


, 


/ i 


V (j/ 2 +4) (^ 2 


^ 1) 



CSch X : 

cschx 



(A16) 



[(1 - v 2 + 2coth 2 x) ± 4wcothx] 



with 



cschx = v/{t+ 1/2) . (A17) 
Where there are zero entries in these equations, a term 
of order Sj Ho/q exists and may be taken into account 
through integration by parts. 

Note that the real part is the contribution to the "im- 
parity of the tensor field and the imaginary part is the 
"i?" -parity contribution. The flat limit is recovered as 
v — > oo, x 0, and v\ — ► kD/Ho. We have numerically 
verified that the underlying expressi on for f alD <5>^ is good 
to ~ 1% for £ > 2 using the code of [Kosowsky] (1998). 



14 



Limber Formulation 

The weak coupling approximation is a gener alization 
of the familia r Limber equation as noted by Jaffe & 
Kamionkowski (1998). In the weak coupling approxima- 



tion, wo are loft with a power spectru m de fined as the in 



tegrall over wavenumher q in equation (|A3|) Tf we c hange 



variables to distance D using the projection relation ( A12) 
we obtain 

Cf x ^tt 2 J dDD A F e J2 Wif | 2 A| (jm) , (A18) 



for a source with a single type of angular dependence S 



(m) 



and logarithmic power spectrum A s 
in I is 



2 (jm) 



The weighting 



Fe=- 



1 



2(^ + 1/2)2 



r[(* + i)/2] 
r[(i + 2)/2] 



tanh( 



25 f + 1/2, 
16 ^D^' 



1 

1* 



(A19) 



The approximation in the second line is for £ 3> 1 where 
contributing modes are well below the curvature scale for 
reasonable cosmologies (v ^> 1). In this same limit, the 
weights also simplify, 



< 0) = o, 



(o) _ 1 



^02 . 

^=±1 



w$ = o, 



w. 



(2) _ 



±22 

(A20) 

Notice that for the spin-2 case, m = 0, ±2 produces a pure 
"E" -parity field and m = ±1 produces a pure "B" par- 
ity field, in contrast to the tight-coupling (delta function 
source approximation) where the power r atios between 
B and E are 0,6,8/13 for m = 0,±1,±2 ( |Hu fc White 



1997). For a simple source S with no angular dependence 
(j = 0, m = 0) and a scalar field (s = 0) on the sky, 
equation ( A18j ) reduces to 
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which is the Fourier Limber equation as derived by Kaiser 
(1992; 1998). 
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